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Abstract. Recently it has been discovered that certain two-qutrit entangled states 
interacting with global and/or multi-local dccohcrcncc undergo distillability sudden 
death (DSD). We investigate this phenomenon for qutrit-qutrit systems interacting 
with statistically independent zero-temperature reservoirs. We show that certain 
initially prepared free-entangled states become bound-entangled in a finite time due to 
the action of Markovian dissipative environment. Moreover, in contrast with local 
dephasing, simple local unitary transformations can completely avoid distillability 
sudden death under amplitude damping. 



PACS numbers: 03.65.Yz, 03.67.Mn, 03.65.Ud, 03.67.Pp 



Submitted to: J. Phys. B: At. Mol. Phys. 



Distillability sudden death in qutrit-qutrit systems under amplitude damping 2 
1. Introduction 

Quantum entanglement is a key resource for quantum information sciences [JJ. As 
the quantum states interact with their environment therefore it is very important to 
understand how an initial entanglement is affected by this interaction. The dynamical 
aspects of entanglement are essential not only for the utilization of entangled states in 
quantum technology but for a deep understanding of the transition from quantum to 
classical world as well. Recently, Yu-Eberly [21 [3] investigated the dynamics of two- 
qubit entanglement under various types of decoherence. They found that although the 
coherence of a single qubit may take infinite time to vanish, nevertheless entanglement 
can disappear at a finite time, a phenomenon named as "entanglement sudden death" 
(ESD). Clearly, such finite-time disentanglement can seriously affect applications of 
entangled states in quantum information processing. The initial report of ESD in two- 
qubit entangled states in a specific model was later explored in wider contexts and 
in higher dimensions of Hilbert space where the qubits are replaced by qutrits (three 
dimensional systems) or qudits (d dimensional systems) [U El [7J [8] . The experimental 
evidences for this effect have been reported for optical setups [9] and atomic ensembles 
[TO] . Recently, it has been observed that [UJ if there is entanglement sudden death in the 
principal system then there is always entanglement sudden birth (ESB) in the reservoirs. 
Interestingly, such entanglement sudden birth may appear earlier, simultaneously or 
later than entanglement sudden death. 

Entanglement of qubit-qubit and qubit-qutrit systems has been completely 
characterized by Peres- Horodecki criterion [T2l [13] . This criterion would imply that 
if a quantum state in Hilbert space of qubit-qubit and qubit-qutrit systems has positive 
partial transpose (PPT), then it is separable. For higher dimensional bipartite systems 
such characterization is not an easy task. Bipartite entangled states are divided 
into free-entangled states and bound-entangled states [HI [15]. Free-entangled states 
can be distilled under local operations and classical communication (LOCC) whereas 
bound-entangled states can not be distilled to pure-state entanglement no matter how 
many copies are available. It was observed that bound entanglement may activate 
teleportation fidelity [16] and manifests the irreversibility in asymptotic manipulation of 
entanglement p2]. Bound entanglement constructed by purely mathematical arguments 
have existence in physical processes [T8] . 

Recently, it was proposed that free-entangled states may be converted into bound- 
entangled states under either multi-local decoherence [19] and/or under combination of 
collective and local dephasing processes [20]. In particular, it was shown that certain 
free-entangled states of qutrit-qutrit systems become non-distillable in a finite time 
under the influence of classical noise. Such behavior has been named as distillability 
sudden death (DSD) [1JJ]. This discovery raises an interesting question that whether such 
behavior is restricted to a certain type of environment or it is also a generic feature like 
entanglement sudden death. We have studied this phenomenon for certain qutrit-qutrit 
states which undergo the simple process of spontaneous emission. We have analyzed the 
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dissipative dynamics of three-level V-type atomic system where spontaneous emission 
may take place from two excited levels to the ground state and direct transition between 
excited levels is not allowed. We have found that distillability sudden death may occur 
in this type of open systems as well. This suggests that DSD might be a generic feature 
such that free-entangled quantum states dynamically convert to bound-entangled states. 
However, we have noticed that unlike local and global dephasing processes, simple local 
unitary transformations can avoid this phenomenon in our specific model. 

This paper is organized as follows. In Sec. El we discuss the physical model and the 
basic equation of motion along with its solution. In Sec. [3l we discuss briefly the idea 
of distillability sudden death and demonstrate the possibility of distillability sudden 
death under amplitude damping. We show that locally equivalent quantum states do 
not exhibit such behavior. We discuss the physical reason for distillability sudden death 
in a particular family of quantum states. We conclude our work in Sec. HI 



2. Dynamics of qutrit-qutrit system under amplitude damping 

2.1. Three-level atom and quantum interference 

We consider a three-level atom in the V configuration. Let |1), |2) be the two 
nondegenerate excited states with transition frequencies to the ground state |3) 
represented by ui, u 2 and electric dipole moments jli, fa respectively. We assume 
that excited states can decay to the ground state by spontaneous emission and a direct 
transition between excited states is not allowed. The spontaneous decay constants of 
the excited states \j), j = 1, 2 are given as 

It is well known that if the dipole moments of these two transitions are parallel, then an 
indirect coupling between levels |1) and |2) can appear due to interaction with vacuum 
[2T] . This coupling between transitions |1) — > |3) and |2) — > |3) is termed as quantum 
interference with cross damping constant 

7i2 = A V7i 72 , (2) 

where 

\tii\ \M 

The parameter (3i represents the mutual orientation of transition dipole moments such 
that for maximum quantum interference f3\ = 1 and for /3j = it (interference) vanishes. 
Figure Q] shows this typical configuration. 

The time evolution of this system-I is given by the master equation [221 123] 

^ = -i[H,p]+A lP , (4) 
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Figure 1. Three level atom is shown in V configuration as system-I. 

where the damping term is given as 

Ti T2 
Ai p = — (2 cr 3 i p 0-13 - oil P - P o-n) + y (2 o- 32 p 0-23 

T12 

- 0"22 P - P 022) + -y (2 O-31 p 0- 23 - 021 P - /3 O21) 
T12 

+ -y (2 0-32 p 0-13 - 0-12 P - po-12) • (5) 

Here ou is the atomic transition operator taking an atom from level |/) to \k). As we 
are interested in the time evolution of initial states due to the spontaneous emission, we 
can ignore the Hamiltonian part (for details see [22} 123]). 




Figure 2. Three level atom (system-II) is shown in V configuration. Two excited 
states |e) and \u) are allowed to make transitions to ground state \g) with decay rate 
7e and j u , respectively. The transition between excited states is dipole forbidden. 



In atomic spectroscopy, the transition dipole moments are usually perpendicular 
and system-I described by Eq.()5]) is difficult to realize. However, it was shown [22] that 
the effects of quantum interference can be simulated to a large degree by three-level 
system-II as shown in Figure |2] with excited states |e), \u) and the ground state \g). 
The corresponding master equation contains the damping operator, 

An p = — (2 o ge p a eg - a ee p - p a ee ) + — (2 a gu p o ug 

- cr uu p - pa uu ) . (6) 
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The main arguments from Ref. [22J supporting this idea are sketched as follows. We 
consider the system-I with 71 = 72 = 7- By taking symmetric and antisymmetric 
superpositions of the excited states 

|a> = ^(|l) + |2», |a) = -L(|l>-|2», 

and the ground state |3) as a new basis in C 3 , the damping term (jSJ) becomes 

7 7 
Aip = -(1 + /3i)(a 3s pa s3 - a ss p- pa ss ) + -(1 - fa) 

x (cr 3a p(T a3 - a aa p - pa aa ) . (7) 

We note that for almost parallel dipole moments (3i ~ 1, the antisymmetric level is 
metastable. The system-II with damping operator (jHJ) is a generalization of this system 
in which we allow the two excited states |e) and \u) to be nondegenerate with arbitrary 
decay rates to the ground state \g). Moreover, the dipole moments can be perpendicular. 
A measure of quantum interference is defined [22} [23] as 

Ai = • (8) 

This expression shows that maximum interference are expected in system-II when the 
level \u) is metastable. It was demonstrated [22] that system-II simulates quantum 
interference in a three-level atom. Moreover, it is simpler model and can be studied 
analytically. We focus on this model and concentrate on the question that whether 
certain free-entangled states undergo distillability sudden death or not. As the effects of 
quantum interference on disentanglement have been studied [23], we would just take a 
fixed value of interference parameter and study our question of interest in what follows. 

2.2. Time evolution of pair of three-level system 

We consider the system of two three-level atoms A and B both of type (II) such that 
they are separated by a distance which is much larger than wavelength of radiation. 
This condition implies that both atoms are independently interacting with their own 
local reservoirs. This case is simpler to analyze and can be solved exactly. The two 
atoms may be initially prepared in an entangled state. The master equation describing 
the dissipative part of dynamics is given as [23] 
dp 



with 



Ap = f (2<p< - aip- pat) + f (2<p< 
-°tuP~ P °tu) + -w (2 o-g e p 0% - af e p - pa\ 



+ y (2 < p or* - a B au p-p a B au ) , (10) 



and where 



crfj = ay <g> I 3 , o-fj = h ® o-ij , i,j = e, u, g. (11) 
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The spontaneous emission of atoms A and B from their excited states |e) to their ground 
states \g) is described by the spontaneous decay rate 7 e whereas 7„ is the spontaneous 
decay rate of excited states \u) to ground states \g). The atomic transition operator 
a ij = takes an atom from its state \j) to state |z). Let us identify |2), |1), and |0) 
be the first excited state, second excited, and ground state of the qutrit, respectively. 
We choose the basis { 1 2, 2), |2,1), |2,0), |1,2), |1,1), |1,0), |0,2), |0,1), |0,0)}. The 
most general solution of the master equation ([9]) has been provided for matrix elements 
of any initial density matrix in Appendix A. From this general solution, it is obvious 
that for 7^ > 0, all initial states dynamically reduce to ground state of both atoms at 
sufficiently large time (t — > oo), which is a separable state i.e. |0,0). For maximum 
interference (7„ = 0), there may be some asymptotic entangled states [23]. However, we 
are only interested in the dynamics with 7 U > 0. 

3. Distillability sudden death 

First we discuss the characterization of bound-entangled states and then demonstrate 
the existence of distillability sudden death. It was shown [15] that bound-entangled 
states have positive partial transpose (PPT) and therefore non-distillable under LOCC. 
In general, there is not a unique criterion to detect all bound-entangled states. Even for 
qutrit-qutrit system, there are various bound-entangled states and a single criterion is 
not capable to detect all of them (see for example [24] and references therein). However, 
one can use the realignment criterion [25] to detect certain bound-entangled states. The 
realignment of a given density matrix is obtained as (p R )ij t ki = Pikji- For a separable 
state p, realignment criterion implies that Hp^H < 1. For a PPT-state, the positive value 
of the quantity \\p \\ — 1 can prove the bound-entangled state. 

Distillability sudden death corresponds to a situation where free-entangled states 
become non-distillable in a finite time under the influence of decoherence [19J. It has 
been shown explicitly that PPT-states are non-distillable (see Ref. [24J and references 
therein). The quantum states with negative partial transpose (NPT) are regarded as 
distillable. There is a conjecture on the existence of NPT-bound entangled (NPT-non- 
distillable) states [26] , however the conclusive evidence is still missing. Therefore if an 
initial NPT-state becomes PPT-state in a finite time under the influence of dynamical 
process, then due to possibility of bound-entangled states, it suffers a transition from 
being distillable to non-distillable. This finite-time phase transition is called distillability 
sudden death (DSD). 

Let us consider a particular family of quantum states given as 

2 . T . . T . a 5 — a , . 

p a = -\m + ){m + \ + -a + + ^-a. 1 (12) 

where 2 < a < 5. In Eq. ([12]) maximally entangled state = l/y/3(\0, 1) + |1,0) + 
|2,2)) is mixed with the separable states a+ = 1/3(|0, 0) (0, 0| + |1,2)(1,2| + |2,1)(2,1|) 
anda_ = 1/3(|1, 1)(1, 1| + 12, 0) (2, 0| + 10, 2) (0, 2|). It was shown [IS] that p a is separable 
for 2 < a < 3, bound-entangled for 3 < a < 4, and free-entangled for 4 < a < 5. 
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Let us take these states with 4 < a < 5. The time evolution of these states can be 
easily determined from the general solution. In particular, we need to determine the 
time evolution of appropriate measure of entanglement for qutrit-qutrit systems. One 
appropriate measure for free-entangled states (NPT-states) is the negativity [27] , which 
is a measure of a state for having negative partial transpose. For positive value of this 
measure, the state will be NPT (hence entangled), whereas when this measure is zero, 
the state is PPT and we can not conclude the entanglement or separability of the state 
until some other measures/steps reveal the status of the state. The negativity is based 
on the partial transpose of a state p. The negativity N(p) is equal to the absolute 
value of the sum of negative eigenvalues of partial transpose of a state p. Once the 
negativity becomes zero, we can easily study the time evolution of realignment criterion 
to determine the existence of bound-entangled states. However, we should keep in mind 
that realignment criterion can not detect all entangled states. 

The partial transpose of the time evolved density matrix p a (t) can have three 
possible negative eigenvalues. However, the expressions for possible negative eigenvalues 
are quite involved, therefore instead of providing them, we prefer to plot them for specific 
choices of parameters a, 7 U , and j e . The negativity for our example is given as 

N(p a (t)) = max[0, —771] + max[0, -772] + max[0, -773] , (13) 

where r/j are the possible negative eigenvalues. 
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Figure 3. The negativity and realignment criterion is plotted against the decay 
parameter j e t for a — 4.2. We have defined 7„ = 0.5 j e . 

Let us now demonstrate that the time evolved density matrix p a (t) do undergo 
distillability sudden death. This could be most easily demonstrated by checking 
entanglement of p a {t) after it becomes PPT. Let us fix the interference parameter 
(3 such that 7 U = 0.5 7 e throughout this paper which means that (3 = 1/3. The 
relation between negativity and realignment criterion is plotted in Figure [3] against 
decay parameter j e t and a specific choice of the single parameter a. Figure [3] shows 
that under amplitude damping an initial free-entangled state becomes bound-entangled 
at a time jetN ~ 0.1826. The entanglement of PPT-state is verified by the positive 
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value of ||p fl (t)|| — 1 in the range 0.1826 < 7 e t < 0.2426. However, the realignment 
criterion fails to detect the possible entanglement after time jetn ~ 0.2426. Hence we 
have demonstrated that free-entangled states exhibit distillability sudden death under 
amplitude damping. 

Figure H] shows the time evolution of negativity and realignment criterion for an 
initial quantum state p a (Eq. ffl2|) ) with a = 4.5 and (3 = 1/3. For this particular 
case, the negativity becomes zero at jet^ ~ 0.6877, however the realignment criterion 
fails to detect the entangled states after time 7 e ^R ~ 0.299. This example shows the 
sensitivity of realignment criterion on parameter a. We have seen that although the 
initial NPT-states (free-entangled) become PPT after a finite time, however we can 
not conclude their separability/entanglement immediately. The PPT-states might be 
entangled suffering distillability sudden death followed by entanglement sudden death. 

0.25 
0.20 
0.15 

o.io 

0.05 



0.00 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Figure 4. The negativity and realignment criterion is plotted against the decay 
parameter j e t for a = 4.5. We have taken 7„ = 0.5 j e . 

Let us now consider the locally equivalent state a a , which may be obtained from 
Eq. (fl"2l by applying the local unitary operation U = I3 <8> 9, with 6 = |0)(1| + |1)(0| + 
|2)(2|. Then we have 

a a = Up a U^ = 2 - (* + | + ja + + ^ ~a_ . (14) 

As the local unitary transformations can not convert separable states into entangled 
states, therefore the transformed separable states are a + = 1/3 ( 1 0, 1)(0, 1| + |1, 2)(1, 2| + 
|2,0)(2,0|), and = 1/3(|1, 0)(1, 0| + |2,1)(2,1| + |0,2)(0,2|). Similarly maximally 
entangled state is also converted into another maximally entangled state |^+) = 
l/v3(|0, 0) + |1, 1) + 1 2, 2) ) . The local unitary transformations do not effect the trace 
and static entanglement of an initial quantum state, however they may have a profound 
influence on the future trajectory of entanglement [7]. In the current situation the 
underlying dynamics is caused by spontaneous emission. For the range of 4 < a < 5, 
the dominant noisy component in Eq. (TT2|) is a + , which contains two doubly excited 
components |1, 2), and |2, 1), hence one could expect that entanglement and distillability 
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may decay faster due to increase of probability of spontaneous emission. What we 
could achieve with local unitary transformations is to modify this component. Now 
the dominant noisy component in Eq. ( 1T41) is <r+ which contains only one doubly 
excited component |1,2), therefore we expect considerable change in the dynamics of 
entanglement and distillability. Indeed, we find the expected intuitive results discussed 
below. 

The possible negative eigenvalues Ai, A 2 , and A 3 of the partial transpose of time 
evolved density matrix a a (t) are again quite involved, therefore instead of writing their 
expressions, we plot these eigenvalues against the decay parameter j e t. In Figure O 
we have plotted the negative eigenvalues against the decay parameter j e t. We observe 
that two of the negative eigenvalues become positive at some finite times, however one 
eigenvalue remains always negative in the range 4 < a < 5. This implies that the state 
a a remains always NPT-entangled until infinity, where it becomes separable. Hence 
distillability sudden death and entanglement sudden death never occur for a a under 
amplitude damping. In contrast the same state exhibits distillability sudden death 
and entanglement sudden death under multi- local dephasing [191 120] • Hence we have 
seen that simple local unitary transformations can completely avoid distillability sudden 
death and entanglement sudden death under amplitude damping. 
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Figure 5. The possible negative eigenvalues are plotted against the decay parameter 
7 e i for a = 4.2. The plot shows that the state remains NPT at all times except 
7 e i = 00. 




It is obvious from Appendix A that for certain classes of initial states, our dynamical 
process evolves within that set. This simply means that all matrix elements which were 
initially zero remain zero. The examples studied in this paper are of such classes. 
Distillability sudden death has been observed so far (including this paper) for a very 
special class of states defined in Eq. (|T2|) and Eq. (fill) [T91 [20] . These states are special in 
the sense that they already contain separable (2 < a < 3), bound-entangled (3 < a < 4) 
and free-entangled (4 < a < 5) states. Our dynamical process decreases all matrix 
elements of initial states except pgg(t). If we start with an initial free-entangled state, 
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then it may happen that such decrement matches with bound-entangled region of this 
family causing distillability sudden death. Further dynamics bring all bound-entangled 
states to separable states. 

Let us consider another important class of states called isotropic states [28]. These 
states are given as 



where < p < 1. The isotropic states have the property that their PPT region is 
always separable (see Ref. [21] and references therein). Therefore, these states do not 
suffer distillability sudden death under any type of decoherence, however they may suffer 
entanglement sudden death [23]. This might indicate that distillability sudden death 
could be related to some particular classes of states which include bound-entangled as 
well as free-entangled states. 

4. Discussions 

We have studied the entanglement dynamics of bipartite systems under amplitude 
damping. In particular, we have focused the distillability sudden death of qutrit-qutrit 
systems. We have found that certain free-entangled qutrit-qutrit states become bound- 
entangled as a consequence of dissipative dynamical process. The explicit demonstration 
of distillability sudden death under the amplitude damping and phase damping suggests 
its generic nature. We have observed that local unitary transformations can avoid 
this phenomenon under amplitude damping. This situation is in contrast with phase 
damping where local unitary transformations can not avoid distillability sudden death. 
Therefore, we conclude that similar to multi-local dephasing case, qutrit-qutrit states 
interacting with statistically independent local reservoirs also evolve into three main 
types of dynamics: (i) They loose their entanglement only at infinity and hence never 
suffer entanglement sudden death and distillability sudden death, (ii) They loose their 
entanglement at a finite time but never undergo distillability sudden death, (iii) They 
first undergo distillability sudden death and then sudden death of entanglement. 
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Appendix A. 

The most general solution of the master equation (Q for an arbitrary initial state is 
given as 



p p = p\y + )(V + \ + 



(1-p) 



(15) 
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Pu(t) = Pn e 
Pi2{t) = Pne 



-2 Jet 



(3 7e+7n) * 
2 
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